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ABSTRACT 

This paper examines the uniqueness and accuracy of the equation 
which describes the transport of charged particles diffusing in a 
random magnetic field parallel to a relatively large guiding field. 

With regard to uniqueness, it is found here that the same coefficient 
of diffusion is obtained by three methods that have apparently led to 
discrepancies in previous work. With regard to accuracy, it is found 
that two corrections must be added to the familiar expression in which 
the diffusive flux is proportional to the gradient of the density. If 
the temporal derivative of the flux is large, the first of these cor- 
rections is important and the diffusion equation must be replaced by 
the telegrapher's equation. If the gradient of the flux is large, the 
second correction is important and the appropriate generalization of the 
diffusion equation is an integral equation analogous to those invoked in 
classical transport theory. Explicit expressions are given for a char- 
acteristic time and a characteristic length which describe the corrections. 



The well known divergence of the coefficient of diffusion, which is 
implied by the quasilinear analysis of pitch angle scattering, does 
not occur if the scattering rate is finite at 90° pitch angle. This 
effect is illustrated by formulas which give the coefficient of diffu- 
sion when the quasilinear expression is perturbed by a variable amount 
of isotropic scattering. 

Subject headings: cosmic rays - hydromagnetics - magnetic fields 
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I. introduction 


The diffusive nature of charged particle propagation in random 
magnetic fields has been recognized for many years. However, there 
are still two unresolved difficulties that obscure the relationship 
between the coefficient of diffusion D and the statistical quantities 
that describe the random fields. In the first place, there is no con- 
cordance on the method of evaluating the Fokker— Planck coefficient 
<Ay 2 >/At which characterizes the scattering of pitch angles relative 
to the average "d.c." field B. The quasilinear approximation (Jokipii 
1966) leads to a simple relationship 


<Ay 2 >/At - (j>{y) = A | v | ^ ^(l-y 2 ) (1) 

in which y is the cosine of the pitch angle and q is the spectral index 
of the power law that gives the mean square amplitude of field fluctua- 
tions at wave number k within an interval dk, Q^Ck^/k^dk, in terms 
of the spectral density Q at a reference wave number k . Here, the 


parameter A, 




can be expressed in terms of the particle rigidity R, velocity V and 
Larmor radius r^ and the spectral parameters Q^, 9 an< ^ sur- 

prising implication of equation (1) that scattering vanishes at y = 0 
for q > 1 has motivated attempts to refine the quasilinear calculation. 
Even the validity of the Fokker-Planck approach has been questioned 
(Lerche 1972; Klimas and Sandri 1971), but a consensus of recent work 
(Jones, Kaiser and Birmingham 1973; Jokipii and Lerche 1973; Volk 1973; 
Owens 1974) agrees that the formalism is applicable but that the coeffi- 
cient <Ay 2 > /At is finite at y » 0. However, the magnitude of this 
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coefficient Is not yet firm. The second unresolved difficulty is that 

different methods of calculating the coefficient of diffusion corresponding 

that the 

to a given Fokker-Planck coefficient give different answers, in spite of the fact/ 
coefficient D is defined uniquely as the ratio of the flux to the gradient 
of the density. This paper is addressed to the second of these diffi- 
culties but not to the first. 

To describe the distribution function for charged particles diffusing 
along a guiding field, three representations have been invoked: (1) A series 
expansion in terms of Legendre polynomials as in classical transport theory 
(Jokipli 19685 Weinberg and Wigner 1958), (2) A series expansion in terms 
of eigenfunctions of the scattering operator (Earl 1973) , (3) A closed 
expression obtained by perturbation techniques (Jokipii 1966; Hasselmann 
and Wibberenz 1970; Kulsrud and Pierce 1969), In previous work, these 
approaches apparently gave different answers. In §11, this conflict is 
resolved by the demonstration that the correct application of either one 
of the first two methods leads to a result identical to that obtained by 
the third method. Thus, the correct expression for the coefficient of 
diffusion is the one derived by Jokipii (1966) . 

This agreement on the description of diffusive streaming motivates 
the consideration , in §111 , of higher order effects which set limits on 
the accuracy of the diffusion equation. These effects lead to the intui- 
tively plausible result that the evolution of the distribution function 
at a point in space and time is controlled not by the local density but 
instead by the average density within a spatial region comparable to the 
mean free path and over a period comparable to the collision time. 



- 5 - 


Jokipii (1971) has reviewed several examples in which coefficients 
of diffusion based upon equation (1) provide a useful description of 
observed phenomena. However, when q - 2, equation (1) implies that D is 
infinite. This divergence is disquieting for two reasons, in the first 
place, for q £ 2 below the point of divergence, the coefficient D is 
very sensitive to variations in the spectral index. Because the index of 
the interplanetary magnetic power spectrum is observed to change with 
time within this region (Sari and Ness 1969; Siscoe et al. 1968), the 
predicted sensitivity implies a greater variability in the local propaga- 
tion of solar and galactic cosmic rays than is observed. In the second 
place, the existence of a physically unreasonable divergence casts doubt 
upon the validity of equation (1) , These doubts stimulated the develop- 
ment of the refined theories mentioned above. The effect of the finite 
scattering at y = 0 predicted by these theories is illustrated, in §111, 
by the explicit calculation of the coefficient of diffusion corresponding 
to a Fokker-Planck coefficient modified » as was suggested by Owens (1974) , by 
the addition of an isotropic scattering term 

$ = A|y | q_1 (l-y 2 ) + H(l-y 2 ), (3) 

where H, the parameter that characterizes the magnitude of the perturba- 
tion, is evidently the coefficient at y = 0 for q > 1. Isotropic scat- 
tering could also arise from a separate effect such as binary Coulomb 
collisions. 

An earlier paper on diffusion (Earl 1973), which introduced the con- 
cept of scattering eigenfunctions and eigenvalues, will be referred to as 
Paper I and equations therein will be designated here by the Roman numeral 
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I. A subsequent paper (Earl 1974), designated hereafter as Paper II, 
treated the coherent mode of particle propagation that replaces diffusion 
when D is infinite. The ideas in Papers I and II will be reviewed below 
and used extensively. 


Particle transport is described by the Boltzmann equation 


9f , 9f _ 1 3 . 3f 
St ^ Sz 29 y ^ 3y 


(4) 


where f {y, z, t) is the particle distribution function, z is distance 

the mean field, and t is time. Because the scattering opera- 
tor on the right hand side of equation (4) has the form specified by 
Sturm-Liouville theory, there is a sequence of eigenfunctions R^Cy) 
and eigenvalues (2 /t^) which satisfy the equation 




(5) 


where is a relaxation time that describes the temporal decay of an 
anisotropy proportional to E R . These functions form an orthogonal set 
in terms of which the distribution function f and the streaming term 
yV(9f/9z) can be expressed as series expansions 

f{y, z, t} = l f R {z, t} ^^{y}, (6) 

pV(3f/3z) = S V j K Of J /9z)c K R K (y}, ( 7 ) 


where the normalizing factors c K defined by equation (II-7) are given 

in table 1 of Paper II. Here, the coefficients f are given by 

K. 

f K = C k£ ] f \ dU * (8) 

and the characteristic velocities V are given by 

JK 

V JK " V KJ " VC J C K J[ ) du * 


(9) 
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When the series representations of the distribution function (eq. [6]) 
and of the streaming term (eq. [7]) are substituted into equation (4), 
the orthogonality of the implies that the coefficients f^ are described 
by a set of coupled differential equations the first four of which are 


_9f 0 

9t 

9fi 


3fj 

'01 97 


= -Vm -r“l -V03 — 3 - 


3f ; 
9z 


I ^ „ \J 3f 0 TT ^ ^ 2 

9t + T 1 ~ Vq1 Tz " Vl2 H * 


9f, 


9z 


9f 2 i f 2 tt j. tt 

_i + _2 = -V l2 -l -V 23 . ... 


3z 


3f 3 , ^3 .. 3^0 tt 3^2 

—a. + -3 = _v 03 — -v 23 — 2 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


9z 9z •** ' 

The matrix formulation of transport theory embodied in these equations is 
entirely equivalent to the Boltzmann equation. 
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II. DIFFUSION REVISITED 

Scattering causes the distribution function to relax toward isotropy. 
An important stage of this evolution is the diffusive regime in which a 
relatively small and nearly constant anisotropy leads to slow temporal 
changes of the isotropic density. In this regime, the distribution func- 
tion f can be approximately represented as the sum of an isotropic com- 
ponent F q = C Q f o which depends upon position and time and two anisotropic 
components which are constant but which depend upon position 

f{]i, z, t} = F q {z, t} + F 1 {p, z} + F 2 {p, z}, (14) 

where F^ is an odd function of p and F 2 is an even function of p. Because 
the form of the odd anisotropy F^ is well known (Jokipii 1966; Kulsrud 
and Pearce 1969; Hasselmann and Wibberenz 1970), this section focuses 
upon relationships among the three methods of analysis listed in the 
introduction. The assumption that the even anisotropy F 2 is negligible 
compared to F q and Fj, which is adopted throughout this section, will 
be confirmed in §111. Nevertheless, it will simplify the treatment that 
follows to require that 

| +1 F 2 {z, p} dp = 0 . (15) 

Evidently, by virtue of equation (8), this condition implies that the 
isotropic term in the series expansion of F 2 is zero. Because Fj^ is an 
odd function, it also satisfies equation (15). 

The odd and even anisotropies are described, respectively, by the 
odd numbered and even numbered coefficients in the eigenfunction expan- 
sion of f (eq. [6]). Thus, if the anisotropies are constant and if the 
even components are negligible compared to the density f Q , then the 
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matrix equations for the odd components take on a form 


f = _ x v 

K K OK 9z 


(16) 


which can be substituted into the gradients appearing on the right hand 
side of equation (10) to yield a diffusion equation for f 


9£o 3 D 9fn 

3t dz 3z 

in which the coefficient of diffusion D can be expressed as 


(17) 


D - xj V 01 2 + t 3 V 03 2 + x 5 V 05 2 + 


(18) 


In Paper I, only the first term of this series was considered. However, 
it is possible to estimate the contribution of the other terms with the 
aid of the identity 

(1/3)V 2 = V 01 2 + V 03 2 + ... (19) 

which results when the series expansion 

p = (/I/VKVoiCiRi + V 03 C 3 R 3 + ...) (20) 

is multiplied by itself and integrated over p. Because t 3 > T 5 > ... , 
the error introduced by retaining only the first term of equation (18) 
is governed by the inequality 

D - T 1 V °1 2 < *3 < 5 V2 - V «l 2 > ‘ ^ 3 (5-2q) (8-3q ) 1 T 1 V 0 1 2 < 2l > 


in which the expression in square brackets represents, for q < 2, an 
upper limit on the fractional error. In the range 1 < q < 2, the maximum 
error of 1.2% predicted by this expression at q = 1.8 is negligible com- 
pared to experimental uncertainties. Only in the unphysical region 
q < -2.1, where the fractional error eventually approaches (-q/18) , is 
the error larger than 10%. 
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According to equation (21) , the error vanishes when q ■= 1 and when 
1 = 2 - the former case, the eigenfunctions are Legendre polynomials 
for which the characteristic velocities V 03 , V 05 , ... are all zero. In 
the latter case, the fractional error is zero because the first term in 
equation (18) approaches infinity while all other terms remain finite. 
These considerations apply even when the Fokker-Planck coefficient is 
described by equation (3), for the eigenfunctions approach Legendre poly- 
nomials when H is large and the first term dominates when H is small. 

Thus, under a variety of circumstances, the flux represented by the higher 
order terms of equation (18) is quantitatively unimportant. However, 
these terms have conceptual significance, because they account for dis- 
crepancies between the analysis given in Paper I and the perturbation 
analysis that follows. 

If equation (14) is substituted into the Boltzmann equation, inde- 
pendent equations are obtained for odd and even components 


at 


3z 


„ 3Fn , „ 3F2 

vV + uV Tz 


2 

i a 


sf 2 

sf 2 
“ at 


(22) 

8F 2 

3y 

3Fj 

at 

• 

(23) 


When equation (22) is Integrated over y from -1 to +1, the scattering 
term on the right hand side contributes nothing because <f> vanishes at 
both limits while the integral over (aF 2 /at) is zero because of equation 


(15). Thus, the integral of the left hand side gives a fundamental re- 
lationship between the time derivative of F q and the gradient of the flux 
S 




= 0 


(24) 



- 11 - 


where 

S 



yFi dy. 


(25) 


Equation (24) holds even if Fj and F 2 are large or varying 
rapidly. Moreover, equations (24) and (25) imply that the temporal 

evolution of F q is controlled by the odd anisotropy alone. The even 
anisotropy influences F^ only indirectly through the effect on Fj implied 
by equations (22) and (23). Because this influence is weak, the odd 
anisotropy will be considered in detail before the even anisotropy is dis- 
cussed. Also, the time derivatives on the right hand sides of equations 
(22) and (23) , which are zero in the diffusive regime, will be neglected 
in this section. However, the role of (3Fi/3t) will be considered in 


Sill. 


If F 2 is negligible compared to F q , equation (23) leads to 

* f 1 ■ -< i ->' 2 > v H° < 26 > 

in which the constant of integration was chosen to satisfy the condition 
<f>{+l} = <(>{-1} = 0. A second integration with the constant chosen to 
yield an odd function gives 

(27) 

The flux 
takes on 

(28) 


_ _ 3F D f p 1-v 2 

F ' ’ - v afj TuT dv 


a result obtained earlier by Hasselmann and Wibberenz (1970). 
S calculated from equation (27) with the aid of equation (25) 
the diffusive form 
3Ff, 


S - -D 


3z 


in which D is given by 
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D . „ dyjT dv = <v’/2)jf' dv - ' < 29 > 

Here, the second equality follows when the order of Integration over y and 
v Is Interchanged. Equation (29) has also been derived by Jokipii (1966), 
Kulsrud and Pearce (1969) and Hasselmann and Wibberenz (1970). 

To link this approach to the method of eigenfunctions discussed above 
and in Papers I and II, expand the odd anisotropy as an eigenfunction 
series 


Fi = Ci + f 3 c 3 R 3 + ... . 

Then, the result of substituting this series into the right hand side of 
equation (23) can be evaluated with the aid of equation (5) while the 
variable y on the left hand side can be replaced by equation (20) to yield 

f§° <Voi Ci Ri + Vq3 c 3 R 3 + ...) = 

(fl/n) Cj Ri + (f 3 /t 3 ) c 3 R 3 + ... (30) 


in which the relationship. F q = has been invoked. The orthogonality 

of eigenfunctions implies that corresponding coefficients on opposite 
sides of equation (30) are equal. Because this equality leads to exactly 
the coefficients specified by equation (16), it follows that the two 
methods give identical predictions for the odd anisotropy. More speci- 
fically, the functional dependence g{y} of Fj upon y 

*1 - - g{v> ff° “ g{y}(S/D) (31) 

can be described either by the integral appearing in equation (27) or 
by an eigenfunction series 

8{y} = V j[ y ^T dV = ^ T 1 V 01 C 1 R 1 + t 3 V 03 C3 r 3 +•••) = [ ^!~) A 3 . (32) 
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Similarly, the coefficient of diffusion expressed in terms of relaxation 
times and characteristic velocities by equation (18) is identical to the 
coefficient given in terms of an integral by equation (29). 

From equations (31) and (32) , it is a simple matter to calculate the 
anisotropy 6^ associated with the component. 


6 1 = 


1 “ 1 4 

max min 

I + I . 
max min 


r ,, 1 3F 0 g(l} Vg{l> S _ r/ . S 

r W„- [4 - ql VF„ • 


(33) 


In these expressions, the parameter g{l) “ [V/(2-q)A] is a characteristic 
length for scattering. The anisotropy is given by the ratio of this length 
to the scale length for spatial variations of the density. The last 
expression, which replaces equation (1-42), reduces, for q = 1, to the 


form 


6=3 (S/VF q ) 

that applies to isotropic scattering. 

From equation (23), it is evident that g{y} is the unique odd func- 
tion that gives a result proportional to y when acted upon by the scat- 
tering operator. Only in the special case that y is an eigenfunction of 
the scattering operator is the odd anisotropy itself proportional to y. 

To make this point explicit, it is useful to expand the distribution 
function in a Legendre series 

f = I a^(z, t} . (34) 

(1-24) , by 

In the diffusive regime, the odd coefficients a are given, as in equation/ 


» . - I&V^JK+II f’ dv Ii-viiid {v) . (35) 

Tt 3z K(K+1) J Q 4>{v} dv K 

In the integration specified by equation (25), the factor y is, of course, 
identical to P x which is orthogonal to all other Legendre polynomials. 
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Consequently, the flux S calculated from equation (34) 


depends only upon the coefficient a^. When a* is evaluated with the aid 
of equation (35), the coef f icient (Va } /3) (^F^ /3z ) _1 is exactly the same as the 
diffusion coefficient that appears after the second equality of equation 
(29). Jokipii (1968) assumed that ai = (3S/V) and that higher order coef- 
ficients are negligible. This assumption erroneously implies that 
Fi = (3S/V)y. Notwithstanding, this form can be substituted into equation 
(26) and the result integrated over y to yield an expression for the dif- 
fusion coefficient 


2V 2 

9/q if>{y} dy 


(36) 


that is correct only for isotropic scattering (Jokipii 1968; also see equation 
[1-37].) In actuality, equation (34) is a valid representation of the 
distribution function whose first term is proportional to y, but the 
higher order terms can not always be neglected. For example, the ratios 
(a 3 /ai) and ( 35 /a j) calculated with the aid of equation ( 35 ) with given 
by equation ( 1 ) are 

13 = Z ^-q) . £5 = 11 (q 2 -l) 

a l 3 ( 6 -q) * 3l 3 ( 8 -q) ( 6 -q) (4-q) * U 7) 

For q = 1, the special case mentioned above applies, and these ratios 
are zero. In general, the higher order Legendre components are very sig- 
nificant. Thus, at q - 2, = -0.58 and ^/aj - 0.23. In contrast, 

equation (32) implies that the higher order scattering eigenfunctions are 
greatly attenuated relative to the first one. Thus , the ratio of scat- 
tering eigenfunctions that corresponds to the first ratio appearing in 
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equation (37) Is 
* 3 V 03 C 3 


T lV()l c l 


<q-l)(2-q)(ll-4q) 

2(7-2q)(5-2q)^ 


(38) 


J_, dy dy 


(40) 


This quantity is zero at q = 1 and q = 2. Between these points, its 
maximum value is 4.8% at q = 1.7. The minor role played in diffusion 
by higher order eigenfunctions is also illustrated by equation (21). 

In principle, the coefficients a^ could also be obtained by solving 
the following set of matrix equations which follow from the orthogonality 
of Legendre polynomials when equation (34) is substituted in the Boltzmann 
equation 

aa K K V Sac-1 | CK+DV 3aK+ 1 . 2 W 1 ,1 v 

at 2k-i at 2k+3 az 4 ;W T „r J 1 ' 

J JK 

where 

•+1 

‘JK —I K dy du J ~\ 

(Compare Weinberg and Wigner [1958, p. 246] .) Streaming is described in 

each of equations (10) - (13) by an infinite number of gradient terms, 

but, in equation (39), the effect of streaming upon the evolution of each 

coefficient depends only upon the gradients of the next higher and the 

next lower coefficients. Offsetting this simplification is the complexity 

of the Legendre representation of the scattering operator displayed in 

matrix form on the right hand side of equation (39). Here, scattering 

has the effect that the temporal derivative of each coefficient is coupled 

not only to itself, as in the eigenfunction representation, but also to 

all other coefficients of the same parity. This comparison brings out the 

important point that, in the matrix representation based upon scattering 

eigenfunctions, the scattering matrix is diagonal. Under most circumstances 
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this property leads to the simplest possible description of the important 
effect of scattering at the cost of slightly increased complexity in the 
description of streaming. In special circumstances, when scattering is 
unimportant compared to streaming, other matrix representations based upon 
linear combinations of eigenfunctions may be appropriate. Because the 
Legendre coefficients in the diffusive regime are given by equation (35) , 
there is no reason to pursue the alternative approach of inverting the 
matrices in equation (39). However, it is important to realize that the 
set of coefficients from which equation (36) was derived does not satisfy 
equation (39) . 

In summary, the perturbation method, the Legendre method, and the 
eigenfunction method give mutually consistent descriptions of the diffu- 
sive odd anisotropy. The same coefficient of diffusion is obtained by 
all three methods. An inappropriate choice of Legendre coefficients that 
has been invoked in previous discussions leads to results that are severely 
limited in applicability. Nevertheless, the Legendre coefficients can 
be chosen to give a representation that is correct but slowly converging. 
The eigenfunction representation converges rapidly, but it offers no ad- 
vantage in the diffusive regime over the simpler perturbation method. 



III. LIMITATIONS OF THE DIFFUSIVE IDEALIZATION 


The diffusion equation is obtained when an approximate expression 
for the flux, equation (28), is substituted into the fundamental equa- 
tion that describes the evolution of the density, equation (2A) . The 
objective of this section is to derive a more comprehensive transport 
equation based upon an improved approximation to the flux. This deriva- 
tion consists of an evaluation of the even anisotropy F 2 as an integral 
of equation (22) followed by a calculation of S which takes into account 
the streaming term yV(9F 2 /9z) and the time derivative (3F^/3t) which 
appear in equation (23) but which were neglected in the derivation of 
equation (28). The streaming effect is fairly straightforward to evaluate 
but the temporal effect is difficult to treat exactly because the depen- 
dence on y of (9Fj/9t) is not necessarily the same as that of Fj itself. 
Because this difficulty is severe only when the anisotropy changes rapidly 
|ti(9Fi/3t)|^|F 1 |, it will be assumed that the same function g{y} describes 
both of these dependences. When this assumption is inapplicable, the 
perturbation method becomes unwieldly, but the method of eigenfunctions 
is efficacious. 

By virtue of equation (24), equation (22) can be recast in the form 


1 9 9F; 

2 9y ^ 9y 


JL 

9z 


( P VF 1 - S ) = - — S 


v 2 r y 1-v 2 




(41) 


is 

where the derivative (SF 2 /9t) has been neglected and where the second 
equality follows from equations (31) and (32), Evidently, F 2 is the 
even function of y that gives a result proportional to yF^ when acted 
upon by the scattering operator. The result of a first integration be- 


tween limits chosen to give the odd function that is required if F 2 is 
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( 42 ) 


to be even, 

vanishes at p * ±1. When the derivative (9F 2 /9y) given by equation (42) 
is integrated over y, the resulting expression for the even anisotropy 
can be written in a form analogous to equation (31) 


F 2 - - h{,> H - h <M> 


(43) 


where the function h{y) is given by 

, 0 A» v dv V} f w dv r 

{y} " Jn *<v) D Jo WTJo 


(v 2 -q 2 ) (l-o 2 ) 
4>{cr } 


da - C 


(44) 


in which the constant of integration C is chosen in such a way that h{y} 

satisfies equation (15) . If <f> is described by equation (1) , 
v 2-q - v 5-2q 


h{y } 


ir 

A J 0 


T^“ 


dv - C 


(45) 


where C can be expressed in terms of the digamma function i|» (Abromowitz 
and Stegun 1964, p. 267), 

C m jJ* V | ~Y~ + dv = X + 3— q} - ~ » (46) 

Equation (43) expresses the important result that the even anisotropy is 
proportional to the gradient of the flux. With the aid of equation (24), 
this relationship can also be put in terms of the temporal derivative of 
the density. 

An appropriate measure of the magnitude of F 2 relative to F^ is the 
ratio of anisotropies 6 2 /^l calculated as in equation (33) for each com- 
ponent separately 
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i* 

«1 


(h{i) - mo}) || 
Ho 

dz 


(h{l} - h{0» 


iljc 

at 


g{i} 


g{i} 


3F f 

dz 


m D(h{l) - h{0}) 1_ a [ >{3-q> - »{(3-q)/2) V j 1_ 

. j 


g{l> 


(4-q) 


( 47 ) 


where L_ = S/(3S/3z). Here, the last expression states that the ratio 

l ) 

of anisotropies is the ratio of a characteristic length for scattering 

to the scale length for spatial variations of the flux . Because Fj 

is smaller than F by a comparable ratio (See eq. [33].), this Implies, 
o 

under most conditions, that F 2 is very small compared to F q . However, 
in the coherent limit D -*■ 00 , F q and Fj are comparable. Under these 
circumstances, F 2 plays an important role that was discussed in Paper 
II. 


To calculate Fj, substitute equation (43) into equation (23) and 


assume 


9Fi , , i as 
at gtu} d at 


to obtain after two integrations 


*i--jT f^(ia-v 2 )v||a + ifJJ g (p} dp - v 0JT' ph(p) dp) 


The flux S calculated from this expression satisfies 


a 2 s 


_l2 + x ii +S=: _ I ) 

x to? at b az 


(49) 


where t and X are, respectively, 

a characteristic time and a characteristic length that will be specified 
below. Equation (49) and 


. (48) 


3F 0 

at 


as 

3z 


( 50 ) 
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which is equation (24) restated here for emphasis, are coupled partial 


differential equations in F q and S that take the place of the diffusion 
equation. Evidently, the corrections for the streaming and temporal change 
of anisotropies are embodied in the first two terms in equation (49) , 
for it reduces to the diffusive form (eq. [28]) when they can be neglected. 

The characteristic time t obtained by performing on the second term 
of equation (48) the integration over p specified by equation (25) is 
(1-u 2 


V 2 f’ (1-v 2 ) dv f 1 , , , V 2 f' , /fPl-v 2 . ' 

- D J 0 J, * {p} dp ■ D J 0 ThT *>, 


,? = (51) 

/ A(5-2q) (2-q) ^ L) 


where the first equality results from the same interchange in the order 
of integration over p and v that was invoked in equation (29) and where 
the second equality involves a further change in the order of integration 
over v and p together with the specification of g{p} by equation (32). 
Similarly, the characteristic length X is given by, 

= £2 (2-q)(4-q) [ 2 } - ip (%(4-q) } - tJ> {ig(10-3q)}j . (52) 

The two new transport parameters, x and X 2 , exhibit the same dependence 
on (2-q)” 1 that leads to the divergence of D at q = 2. (See eq. [29].) 

The first correction term in equation (49) can be evaluated in terms 
of F q with the aid of equation (50) while the second term can be evaluated 
with the aid of equation (28) provided that the density changes slowly. 

Thus, equation (49) reduces to 
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3z3t 


(53) 


in which the relative magnitudes of the corrections for the even aniso- 
tropy and temporal changes are respectively proportional to X 2 and tD. 
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In Figure 1, where these parameters are plotted as a function of q, the 
quantity 


tD = 


(V 2 /A) 


(54) 


(5-2q) (2-q) 2 * 

which is 3.75 times larger than X 2 at q - 1, becomes infinitely larger 
at q = 2. Because of this divergence, the temporal effect dominates 
over the streaming correction throughout the range 1 < q < 2. Thus, it 
is appropriate to consider the limit X 2 << tD in which equations (49) 
and (50) reduce to the telegrapher's equation 
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where the parameter 


(D/t) 


y ^-2q) 

(4-q) 


(56) 


is a characteristic velocity for the coherent propagation of density 
inhomogeneities. This parameter is the same as the velocity V 01 (See 
table II-l.) that appeared in the telegrapher's equation obtained in 
Papers I and II. Moreover, It was demonstrated in §11 that the coeffi- 
cient of diffusion appearing here is only slightly different from the one 
obtained in Paper I. Thus, if the term proportional to X 2 is neglected, 
the implications of equations (49) and (50) are qualitatively identical 
and quantitatively similar to those brought out in Paper II - §111 under 
the assumption that only the components f and ^ need be considered in 
the eigenfunction expansion of f. 

To make clear the significance of X, equation (49) can be solved 
for S to yield an expression. 



in which the Green's function within the brackets describes the effect 
at (z, t) of a localized, impulsive source of flux at (z o> t Q ) and in 
which the quantity D(3F /3 z) , evaluated at (z ,t ), plays the role of a 
source function. Because an integral equation for F q is obtained when 
this expression for S is substituted in equation (50) , this reformulation 
does not simplify the problem. However, equation (57) leads to the con- 
ceptually significant implication that S is influenced not by local con- 
ditions but instead by conditions over a finite temporal interval char- 
acterized by t and over a finite spatial interval characterized by X. 

Thus , the present approach embodies a non-local quality of the transport 
equations similar to the attribute emphasized by Klimas and Sandri (1973) . 
There is a temptation to identify the spatial dispersion implied by this 
characteristic with the dispersive spreading of coherent pulses that was 
discussed in Paper II - §IV. However, the latter phenomenon includes a 
dispersive effect which arises because the dependence of F^ on y differs 
from that of (aF^at). In Paper II, this difference was taken into 
account approximately by allowing the evolution of the third odd compo- 
nent f 3 to be different from that of the first odd component fj. Be- 
cause equation (57) takes into account neither this effect nor temporal 
variations of F 2 , the method of eigenfunctions is to be preferred for 
the analysis of problems which involve rapid temporal variations. 

On the other hand, a useful description of the steady state is ob- 
tained when the integration over t Q in equation (57) is carried out with 
the source function assumed to be independent of t Q , 



This expression, with S held constant by virtue of equation (50), is an 
integral equation for that describes more accurately than does the 
diffusion equation the steady state spatial density profile in the vicinity 
of localized sources or of abrupt changes in D. The analysis of such 
situations is a major concern of classical transport theory, but this 
objective is less important in the present context than the analysis of 
transients. 

The diffusion equation is an approximation that applies only when 
temporal and spatial variations of the flux are not pronounced. If sig- 
nificant changes occur within one characteristic time t, then the tele- 
grapher's equation (eq. [55]) describes a generalization of diffusion 
that incorporates localized disturbances propagating coherently but with- 
out dispersion. The dispersive spreading of these disturbances involves 
complex behavior of the anisotropy that is best described in terms 

of eigenfunctions. When temporal changes are gradual but significant 
variations occur over a distance X, then the appropriate generalization 
of the diffusion equation Is equation (58) treated as an integral equa- 


tion for the density. 
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IV. THE COEFFICIENT OF DIFFUSION 

Of the three transport parameters, the coefficient D is the most 
important because it appears by itself in the diffusive approximation to 
the flux. The parameters x and X appear only in the correction terms of 
a higher order approximation, but even here D 

plays an important role. Thus, equation (29) takes on special signifi- 
cance as the unique relationship between D and the Fokker-Planck coeffi- 
cient <p . Because this relationship involves only a single integration 
over y, the coefficient D is readily evaluated provided that the integral 
exists . 

The divergence of the integral when <p is given by equation (1) is 
crucially related to the behavior of the Fokker Planck coefficient near 
y = 0. Specifically, when both the coefficient and its first derivative 
vanish at y = 0, then D is infinite. These conditions occur when q > 2. 
F° r 1 < 9 < 2, <J) has a cusp at y = 0 in which the coefficient vanishes, 
but D is finite nevertheless. For q < 1, D is well defined even 
though <j> is infinite at y = 0. 

The coefficient D is always finite if p is given by equation (3) . 
Thus, the qualitative effect of scattering at y = 0 is to eliminate the 
divergence of D. In this section, the detailed nature of this effect is 
illustrated by analytic expressions which give D as a function of the 
parameters g = (H/A) and q. 

Numerical integration of equation (29) is the most generally applic- 
able method of computing D. In fact, it is the only practical method 
when the dependence of 0 upon y is very complicated as in the results of 
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Volk, et al (1974). Thus, the expressions presented here, although 
they apply to a specialized form of 4>, may be of value as an easily im- 
plemented test of computer routines. 

For q > 1, the parameter H is the Fokker Planck coefficient at y = 0. 
The theories of Jones, Kaiser and Birmingham (1973) and of Volk (1974) 
predict the numerical magnitude of this quantity, but the former also 
predicts a dependence of <Ay 2 >/At upon y that can be specified in terms 
of Bessel functions while the latter approximates the dependence by in- 
troducing a critical value of y below which the coefficient is constant 
and above which the quasilinear result applies. Both of these dependences 
are described with sufficient precision by equation (3) when H << A. 

However the inaccuracies that occur when this condition is not satisfied 
are of minor significance compared to those resulting from the neglect 
of scattering at y = 0. Consequently, to provide insight, 

the analysis will be carried out without regard to this limitation. More- 
over, the form of equation (3) is convenient because the eigenfunctions 
reduce to Legendre polynomials when H >> A. In this limit, the quasi- 
linear term can be regarded as a small perturbation of isotropic scattering. 
Note that the divergence in D would also be removed by a perturbation 
with zero value but finite slope at the origin. 


When equation (3) is substituted for <j>, equation (29) becomes 

„ _ v 2 if 1 i-v 2 , __ v 2 i f _ x <^-q)/(q- 1 ) 

° = A 2J o g + v q-l dV = A 2(q-l)J 0 x+8 


dx 


(59) 
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where x = \ Evidently, D is given by a dimensionless function of B 

multiplied by the parameter (V 2 /A) specified by equation (11-86). Although 
a general expression for the integral over x can be written in terms of 
hypergeometric functions (Gradshteyn and Ryzhik 1965, p. 285, eq. [3.194-5]), 
it is more appropriate in the present context to consider the expressions 
that apply when q takes on certain specific values for which the integral 
reduces to elementary functions. These expressions are given in table 1 
for five values of q. In the limit 3 >> 1 where isotropic scattering is 
dominant, the following series expansion In B' 1 , obtained from equation 
(59) by representing the denominator as a geometric series, is applicable: 

DA _ 1_ 1 1 , 1 1 . (60) 

W 33 " q(q+2) Q 7 + (4q 2 -l) B 3 * * * 


Here, the first term embodies the limit 


D 


vi 

3H 


(61) 


that applies in the case of pure isotropic scattering. Limiting forms 
for B ■+ 0 are readily obtained for specific values of q, but they are 
not significantly less complicated than the formulas in table 1. 

In Figure 2, the dimensionless quantity (DA/V 2 ) is plotted as a 

(also considered by Owens [1974]) 

function of 8. Here, the curves for q = 1.5 fend q = 1.75 approach, in 
the limit B -*■ 0, the finite values of (DA/V 2 ) predicted by equation (29) 
which are 0.80 and 1.78, respectively. In contrast, the curves for q = 2 
and q = 3 approach infinity as M 0. Thus, in this limit, the coefficient 
of diffusion depends, for q i 2, upon both H and A whereas it depends, 
for q < 2, only upon A. All of the curves in Figure 1 converge in the 
limit 6 >> 1 to the dashed line specified by equation (61) . In particular, 
this line represents a limit in which, for q >> 2, D does not depend upon 
A. In this limit where the magnetic power spectrum is very steep, the 
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quasilinear term plays a minor role and the formula for isotropic scattering 
applies. In any case, equation (61) gives an absolute upper limit on the 
value of D. 

Direct observations of interplanetary magnetic power spectra (Siscoe 
et al. 1968; Sari and Ness 1969) , suggest spectral indices in the range 
1.5 < q < 2. and give fairly accurate values for the parameter A. Unfor- 
tunately, while existing theories agree that scattering occurs at y - 0, 
they disagree sharply on the predicted magnitude of 8. On the one hand, 
Jones, Kaiser and Birmingham (1973) give 

8c^<6B2> 0 2 (62) 

where a = 1 and where 6B 2 is the mean square amplitude of perpendicular 
fluctuations of the field. This result puts (DA/V 2 ) in the 

center of Figure 2 where the coefficient of diffusion is sensitively de- 
pendent upon both 8 and A. On the other hand, it follows from the work 
of Volk (1973, eq. [52]) that 
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-(=) 
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= r— 22— i 

L 8(3q-l)J 


2/(3-q) 


2tt ( q -1) <AB Z > 


2-> T (q-l)/(3-q) 


_ <AB 2 > 

= .95 — 


B' 


(63) 


where <AB 2 > is the mean square amplitude of fluctuations above the reso- 
nant wavenumber k reg = (r^) -1 and where the second expression involving 
a numerical coefficient applies to the same case, q * 2, considered by 
Jones et al. These two expressions have very different implications. 

In the first place, their coefficients and exponents are such that for 
a given fractional amplitude, which would typically be less than one, 
the value of £ predicted by the latter expression is smaller than that 
predicted by the former one. In the second place, the fractional amplitude 
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dpp BciiTilng m equation (63) is comparatively small because it represents 
an integration over a limited range of wavenumbers. Because of these 
considerations, the formulation of Volk gives smaller values of B and 
larger coefficients of diffusion than the calculation of Jones et al< 

* In particular, the coherent propagation of kilovolt solar elec*" 
trons , whose scattering is controlled by a relatively steep portion of 
the magnetic power spectrum and whose resonant wavenumbers are large, 
could occur, under favorable conditions, according to equation (63) but 
is virtually ruled out by equation (62). (See Paper II for a discussion 
of this coherent mode which may occur when D is very large. Lin [1974] 
has summarized the observations of coherent effects.) On the other hand, 
in the example considered by Volk (1973), the propagation of 100 MeV 
protons is essentially diffusive. Evidently the discrepancies in the 
evaluation of 6 must be resolved before these issues can be clarified. 

This research was supported by the National Aeronautics and Space 
Administration under grant NGR 21-002-066. 
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FIGURE CAPTIONS 

Figure 1 . The correction to the flux arising from temporal changes 
of the odd anisotropy is proportional to tD. The correction arising 
from gradients of the even anisotropy is proportional to X 2 . For 
q > - 0.7, the first of these corrections is larger than the second. 
Figure 2. As the spectral index q is increased, the dimensionless 
parameter (DA/V 2 ) becomes more sensitively dependent upon g = (H/A). 



TABLE 1 


FORMULAS FOR THE COEFFICIENT OF DIFFUSION 
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